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Abstract. Collective Shapiro Phase Shift on the period of the pulsars due to a dark 
object (a passing MACHO or a planet) along the line-of-sight of one or more pulsar is a 
, formidable gravitational tool to discover dark matter. 

oo : 

Q\ • We already noted that the presence of a few negative pulsar periods, possibly due to such 
' a Shapiro delay, is roughly consistent with independent estimates of MACHOs observed 
^ by microlensing in our Galaxy. Here we update our study and suggest to verify and to 
f*i , calibrate the Shapiro phase delay to observe on ecliptic main known planetary and solar 
C3 gravitational fields. Once the test is probed we propose to use, in a collective way, a 
^ ■ subsample of pulsars on the ecliptic plane to monitor our solar system for discovering 
' any heavy unknown object, either dark MACHOs or unknown planets. The importance 
5r i °f such unknown possible planetary components and their eventual near encounter with 
Earth in the far past (or future) should be not neglected : the same life evolution (or 
extinction) might be marked by such a rare event. 
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1. Introduction 



Gravity deflects space-time. For this reason matter and energy, while bending geodesies, 
force other masses to move along Keplerian trajectories. Indeed gravity, by equivalence 
principle and better by General Relativity, bends also massless photons. Deflection of 
star's lights by Sun, nearly eighty years ago, overthrown newtonian relativity leading 
to the most popular triumph of Einstein's relativity. It is the gravitational bending of 



the light that gave life to gravitational microlensing (Paczyhski 1986 1991 ; Griest 1991 ; 
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Alcock 1995), one of the best ways to see "matter in the dark". However gravity shapes 
also time rate. Indeed clocks are slowed in presence of heavy bodies (Earth, etc.). The 
well known gravitational redshift is (as well as light bending) a cornerstone of General 
Relativity. Time is slowed down near pulsars and totally frozen when observed, far away, 
near the Schwarzschild radius of a black hole. 

The variable gravitational redshift on a fixed period due to the source motion with 



respect to the gravitating body is just "Shapiro Phase Shift" (Shapiro 1964 ). It has 
been predicted on early 1964 and observed, by delayed radar echoes grazing the Sun and 
reflected by Mercury (1967) and Venus (1970), by Shapiro himself. We have not (yet) 
a powerful radar to inspect far away dark planets or MACHOs in our galaxy. However 
nature does offer very precise clocks spread far in our galaxy: pulsars. The period of 
pulsars is stable, million or billion times better than our commercial clocks, competitive 
with the best laboratory ones on Earth. The stability of the timing of pulsars is granted 
by the huge inertial momentum of the pulsar, its consequent large angular momentum 
and the negligible angular momentum losses. Therefore pulsars are the timing candles 
in our Shapiro delay search: their beating slowdown and speed up may inform us on the 
dark objects crossing along the line-of-sight. 

While geodesic accelerations suffer the gravitational field gradients (proportional to the 
inverse of square distance) at a given point, Shapiro Phase Shift (SPS) is a path integral of 
the field (which decreases as the inverse of the distances): SPS "records" the gravitational 
fields from the source to the observer and the delay is a cumulative effect. Indeed there 
is both a geometrical delay (which is proportional to gravitational fields) as well as a 
gravitational Shapiro delay which is roughly the path integral of the field. The consequent 
Gravitational Shapiro delay is a logarithmic function of the impact parameter and it is 
usually the dominant one at large impact parameter values of body encounters. This 
effect increases, nearly by a factor of one hundred, the characteristic impact distance of 
a detectable phase delay with respect to a microlensing event of magnification of the 
brightness of stars. Therefore, even the few hundreds pulsars in our galaxy offer a large 



enough sample to measure a few SPS delay a year (Fargion & Conversano 1996 , 1997 ). 
Their imprint is a long-term "anomalous" negative contribute on the time derivative of 
the period of pulsars, superimposed to the slow-down positive one. Depending on the 
strength of the event, the total P can be still positive but lower. Moreover a sub-sample 
of pulsars close to the ecliptic projection may monitor the Solar System planets (Jupiter, 
Mars ...) as well as the solar SPS, allowing a quantitative calibration and verification of 
the SPS delay. The same technique may be deploied to observe far hypothetical unknown 
planets which may pollute at the far periphery of the Solar System. Their eventual rare 
infall and near encounter with the Earth in the past might have left important geological 
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imprints (Fargion & Dar Fargion & Dar 1998); strong tidal forces might induce Tsunami 
and volcanic activities as well as life extinction. Shapiro Phase Delay may also reveal 
MACHOs in our Solar System neighbors while a collective correlation of the period 
derivatives may track the secret trajectories of far dark bodies in the Space. 



2. Review of the Shapiro effect on pulsars 



The Shapiro Phase Shift on Pulsars (SPSP) is the Shapiro effect (Shapiro 1964 ) which 
influences the time of arrivals of the signal of the pulsars. 

The application of the SPSP to the search of dark matter was first proposed by J. 



Schneider ( 1990 ). Limiting to the case of strong lensing events, < b >~ Re, where 
< b > is the mean impact parameter estimate based on visible matter, he constrained 
the observation of the phenomenon to the signals of extragalactic pulsars. Larchenkova & 
Doroshenko ( 1995| ) probably detected an event of SPSP towards the galactic anticenter 
region. The authors claim the existence of a massive black hole of ~ 330 Mq producing a 
strong lensing event. Wex, Gil & Sendyk ( 1996| ) investigated the possibility of identifying 
MACHOs in the galactic center region by means of SPSP both for strong and weak 
lensing events. According to the authors the probability for observing an SPSP with lens 
mass between 1 and 20 Mq is quite reasonable but, for weak lensing events, it should not 
be possible to extract any information on the parameters of the lens (mass, transverse 
velocity ...) without making assumptions about both the intrinsic frequency v and the 
intrinsic frequency derivative v. 



We also analyzed the SPSP effect in our previous works (Fargion & Conversano |1996] , 



1997). SPSP may be responsible of the time derivative deviation and may explain the 



presence of a few negative pulse derivatives among 706 known pulsars. We estimated 
the mean optical depth for SPSP events for weak lensing (high impact parameters) on 
the basis of present models of the distribution of luminous matter in our galaxy. We 
found that the SPSP may play a role in P deviations for just a ~ 0.7 ev yr~ l event 
rate at a |P| < 6.9 • 10~ 14 s.r' level for massive objects of M ~ Mq. We expect a 
higher event rate for masses M ~ 10 -2 Mq (~ 7 evyr^ 1 ) but at lower intensities of 
\P\ < 6.9 • 10~ 15 s s _1 . Events toward the galactic anticenter should be rarer to observe, 
as a result of the rapid decrease of the density of galactic matter with the distance from 
the center. Moreover, the higher the typical mass of massive objects the lower the event 
rate for SPSP should be expected. On the basis of these simple probability estimates 
we believe the event detected by Larchenkova & Doroshenko be rather a weak lensing 
event produced by a MACHO of mass much smaller than ~ 330 Mq. We also showed 
how analogous collective SPSP in globular clusters would also be possible. 
The fundamental equations for SPSP can be summarized as follows. A massive object 
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passing by the line-of-sight of a pulsar, with transverse velocity v± at a distance Dd from 
the observer (see figure |l|), produces a Shapiro time delay. This is made of two terms: 
a first geometrical one, At geo = 5 s - ( Vu 2 + 4 — u ) , which is the result of the angular 
deflection of the path and it is rapidly vanishing at large impact values: At geo | u;s>1 = ^rjfr- 
The second term is the gravitational redshift due to the deflector field, ~ J — , all 
along the wave trajectory and scales as the logarithm of the impact parameter u: 
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where D s is the distance from the observer to the source and u — b / Re is the 
dimcnsionless impact parameter of the MACHO from the line-of-sight of the source, in 
units of the Einstein radius 
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where u m in is the minimum impact parameter of the MACHO (see figure |J). The as- 
sumption of 10 2 as a characteristic value for u rn i n offers a square geometrical probability 
amplification (wu 2 ) to reveal a massive object, making the SPSP on 700 pulsars statisti- 
cally comparable with the ordinary microlensing of intensity magnification on a sample 
of 7 millions of stars. 

The outstanding signature on the period of a pulsar should be a negative contribute to 
the total time derivative during the second half of the SPSP event. The evolution of the 
P deviation is: 
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gives the strength of the event and 

(t) 
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its time evolution. The function F(t) is centered at the instant t = of closest approach. 
It has a maximum and a symmetric negative minimum for the values =p t c (as shown in 



D. Fargion & R. Conversano: Title Suppressed Due to Excessive Length 5 
figure U for values u m i n = 1 and u m i n > 10 ). 

Competitive gravitational delays in residual periods of the pulsars may occur due to 
Doppler shifts induced by gravitational deflection of both unbounded objects (planets, 
stars, black holes, ...) and binary bounded objects near the pulsar (the former being 
a step-like residual signal in the period of the pulsar while the latter being a periodic 
sinousoidal one). In this context the SPSP leads to a recognizable, episodic, timely 
bell-shaped period variation. SPSP period derivative should be added to the average 
positive derivative (10 -14 ss" 1 ) related to the intrinsic pulsar spin-down. Our first 



inspection (Fargion & Conversano 199£, 1997 ) into the catalogue of 558 pulsars (Taylor 
et al. 1993 ) found out some "anomalous" negative period derivatives. In table [l] we 
show the comparison between the two editions of the catalogues of the known pulsars 



(Taylor et al. 1993 , 1995 ), with the epoch of each measurement. We note that the pulsar 
B0021-72D, which had no measured value in the older catalogue, has a negative period 
derivative with low absolute value. The pulsar B0021-72C, which resides in the same 
globular cluster, has different values of P even if the epoch of measurement is the same 
(notwithstanding this still negative). We believe the newer value of P be consistent, 
being similar to the value of B0021-72D and, maybe, being evidence of a collective SPSP 
within the globular cluster. Particular attention must be dedicated to the pulsar B1813- 
26. It does not belong to any globular cluster nor to any binary system but increases its 
negative P = (-3 ± 3) 1CT 16 s s" 1 to a positive P = (+6.65 ± 0.07) lO" 17 s s^ 1 after a 
time interval of ~ 17.5 yr. This may indicate the second phase of an SPSP event, the 
total P getting again positive values, as well as just an improvement of the precision of 
the measure. 

The reported values of P of the other pulsars in table [j] are equal in both the catalogues. 
On the basis of these data we cannot produce a more precise statistic on the time develop- 
ment of the P, in this sense only an on-line frequently updated database should be useful. 



3. Shapiro Delay in dark planet search close to the Solar System 

Many sources of SPSP can affect the time derivative of the period of the pulsars: lo- 
cal companions in binary systems, encounters of MACHOs along the trajectory of the 



signal (Paczyhski 1986, 1991; Griest 1991; Alcock 1995), gravitational waves, collective 



SPSP in globular clusters (Fargion & Conversano 1997). Moreover any massive object 
in Solar System can, in principle, be a source of SPSP. In this section we investigate on 
the possibility to use the SPSP effect to discover dark planets and any other massive 
object close to our Solar System. In the case of planets bounded to our Solar System we 
expect their orbits to lay within low ecliptic latitudes. This defines a a preferring zone 
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of observation close to the ecliptic. We can choose the array of pulsars laying along the 
ecliptic trajectory within an angular radius of ±6.38° (so we can include the pulsars of 
both the galactic center and anticenter). Figure || and table ^ show this subsample. We 
count 73 pulsars which should be continuously monitored. 

There are two main contributes of P-deviation which affect any observed pulsar: a 
Doppler effect D s / c (as a result of both the relative motion of the pulsar and the 
revolution of the Earth) and a Shapiro effect (as a result of the gravitational action of 
the deflector on the optical path of the signal of the pulsar). 
The total time delay due to the two main contributes is: 

A t = A tuoppler + A t(3rav (7) 

We take our attention to the gravitational one. For nearby deflectors the gravitational 
time delay is given by the three-dimensional vectorial formula: 



A t G rav = — In 
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where we adopted the same notation used in section |^. Dd, Dds and D s = Dd + 
Dds are, respectively, the position vectors of the deflector relative to the observer on 
Earth, the pulsar relative to the deflector and the pulsar relative to the observer; Dd = 
\Dd\, D ds = | -Dds | and D s = \D d + D ds \ are the modules and ip = D s D d (figure §). 
When ip — ► 0, cos-0 ~ — 1 and Dd s — D s — Dd- Consequently, after expanding the 
denominator in equation ||, we find f — cos ip ~ (D s D d ) 2 / 2 ~ (6 + x + ) 2 / (2 Dd 2 ) = 
Re 2 {Vu 2 + 4 + u) 2 /(8 D d 2 ) (figures and §) and equation | tends to equation [1} We 
neglect the contribute of the geometrical deviation due to the fact that we are interested 
mainly in weak lensing (large u). The residual period time derivative can be found by 
time differentiating equation ||. 

As the dark planet search involves deflector masses much lower than stellar ones, the 
resulting P-deviations are much lower than what expected in galactic dark matter search. 
Therefore a comprehensive analysis of what are all the possible sources of SPSP in our 
Solar System that can produce SPSP P-deviations is needed. 

The removal of all of these effects, together with the improvement of the threshold of 
sensitivity, is important for the SPSP technique in dark planet search to have success 
in Solar System. Here we propose two ways to achieve this: («) the extraction of the 
frequency component corresponding to annual modulations on the timing; (ii) the cross- 
correlation of the signals coming from an array of pulsars laying along the projection of 
ecliptic toward the galaxy. 
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3.1. Annual Re modulation 

The low-intensity P-deviation scenario in Solar System is quite complex. The Sun and all 
the known planets produce P-deviations higher then or comparable to those expected in 
dark planet search. Even the revolution of the Moon around the Earth and the rotation 
of the Earth itself produce both Doppler and SPSP P-deviations. The evaluation of the 
time evolution of all of these contributes is important to ensure both a calibration and 
the removal of known additive disturbs in using SPSP for dark planet search. All of these 
P-deviations are periodic. Planets in Solar System produce periodic SPSP P-deviations 
with their own periods of revolution, the Sun an annual one, the Moon a monthly one 
and the rotation of the Earth a daily one. The revolution of the Earth produces also an 
annual Doppler P-deviation (the mutual attraction of known planets in Solar System, 
which perturbs Earth trajectory, should be included in this term), while the rotation of 
the Earth a Doppler daily one. 

The variation of the mutual distances among the planets causes the Einstein radius Re to 
change. If we refer the observations to a geocentric observer, the variation is modulated 
by the annual revolution of Earth. The relative variation of the Einstein radius is: 

ARe 1 AAz 1 ADds (Q) 
Re ~ 2 D d 2 D ds U 

where we do not consider the term AD S / D s , being this negligible. We can distinguish 

three cases: 

1. Dd — — D s . This is the most probable case, where the deflector is located nearly 
half-way from the observer and the source. The Einstein radius takes the general form 
of equation |2[ Relative variations of few astronomical units do not affect significantly 
the signal coming from a pulsar located at distances D s ~ Kpc. In these cases planets 
orbiting in distant stellar systems should not give detectable modulations due to the 
orbit of the planet. 

2. Dds <^ Dd — D s . In this case the deflector is close to the pulsar. The Einstein radius 
depends on the deflector-source distance only, Re — v 2 r$ Dds ■ If the pulsar were 
bound in a binary system its signal is modulated at the frequency of rotation when 
ADds gets comparable with Dds- 

3. Dd -C Dds — D s . This is the case when the deflector is very close to our Solar 
System. The Einstein radius takes the form: Re — \J 2 rg Dd , depending only on the 
deflector distance. If the deflector were a dark planet bounded to our Solar System, 
relative variations of few astronomical units should be significant. The revolution of 
Earth itself should infer to the signal an annual modulation as a result of the cyclic 
variation of Dd- As an example, figure]^ shows the typical geocentric orbit of Jupiter, 
together with the regular orbit of the Sun. The mutual Earth- Jupiter distance, Dd, 
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varies periodically nearly twelve times during a whole revolution of Jupiter. The 
period of this variation is exactly one year. Therefore we expect that the P-deviation 
have two periodic components of modulation: a long-term one (with a period of nearly 
twelve years), as a result of the revolution of Jupiter, and a short-term one (with a 
period of one year), as a result of the revolution of Earth. Any massive object passing 
close to our Solar System should infer to the time arrivals of the pulsar a short-term 
Re modulation at a period of one year and, if bounded to Solar System, a long-term 
Re modulation at its own period of revolution. . 

The last item gives an answer to how a close SPSP event can be recognized from a far 
one. -) First, the annual frequency component must be extracted from the timing of the 
pulsar. From this component both the Doppler effect and the solar SPSP, mainly, must 
be subtracted. Now, the resulting signal is still affected by the short-term i?^-modulatcd 
components of some other known planets in our Solar System. The key point is now to 
recognize which planets. 

-) As the SPSP modulations due to the other known planets can be easily calculated, one 
should check if long-term frequency components are present in the residual timing of the 
pulsar. If a long-term component at a known frequency, say, that of Jupiter, is present, 
the corresponding short-term Re component must also be subtracted from the annual 
component. Therefore, as any known planet produces its own characteristic long-term 
component, one can, in principle, recognize the presence of SPSP contributes from known 
planets and subtract them. 

-) Any residual annual frequency component, free from both Doppler and known-SPSP 
effects, if still present, could reveal a massive nearby object. 

Now, what is the sensitivity? We can compare the change in Shapiro time delay, as a 
result of a change in Dd (SdAt grav ) with the P deviation. The procedure to calculate P 



is (Fargion & Conversano 1997 ): 

At grav ® ^ 

du dt 



p ^dAt^du m 



where 



u(t) = u min \l 1 + I — I , (11) 



while 5dAt„ rav is 



5 d At grav ~ ^- AD d . (12) 

du oD d 



Comparing equations |10| and [12| and approximating Re — \/2rsDdi we find S d At grav as 
a function of P and the relative variation ADd/ Dd- 

6 d At grav = At° ■ F d (t) (13) 
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9 



where 



(14) 



and 



F d {t) = (£) F(t) . (15) 
5. 2. Multiple SPSP on an array of pulsars 

Any massive object either orbiting around or passing by close to our Solar System may 
slow down the signal of more than one pulsar at the same time. In fact, the projection 
of its impact parameter subtends an angle wider than that subtended in the case of a 
far deflector. In this case we can use the whole array of known pulsars both to increase 
the detectability threshold and to trace the trajectory of the object via a collective 
statistics. With what threshold of detectability we can hope to achieve this? We can 
choose the maximum acceptable angular impact parameter as, say, the latitude width 
of the subsample in table || (Ab = 2 ■ 6.38°). With this choice we can calculate what 
P-deviation a planet passing by one edge of the subsample generates on a pulsar laying 
on the other edge. The linear impact parameter REUmin becomes 2D d tan(6.38°) and 
the P-deviation (equation |^) produced by the planet is: 

P T = ^ (16) 

T 2D d tan(6.38°) 1 ' 

We can replace the dependence from (3± with the distance Dd of the deflector from the 
observer: 

0± = J® ~ (17) 



where the distance r of the deflector from the Sun has been approximated to Dd for large 
distances, r© = 1 AU and /3® = 9.96 • 10~ 5 . In this case equation [l6| becomes: 

— - 

Ab r S an in _22 ( M \ ( D d \ 2 -1 



Pt- —, r 3 =8.8-10-^ £ ) I ——7— ] ss- L . (18) 

2r e tan(6.38°) / \10- 10 M Q ) \IAU ) K ' 

\~) 

Equation |l^ can be used to determine the minimum mass of a planet that can perturb the 
signal of a pulsar (with almost the maximum impact parameter allowable within the two 
edges) of the subsample at a Pp-level deviation (see figure ||). Then, the characteristic 
time of the Px-level deviation is: 

= P^W = 2P,tan(6.38°) fDd V g 

Figure ^ shows the mass M of the deflector plotted as a function of its distance Dd for 
many threshold levels Pt- Only planets with mass greater than ~ 1CP 3 M Q can produce 
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the actually-observed average P-deviation, P T ~ 10 ~ 14 ss^ 1 , across the two edges of 
the subsample. To detect mini-planets (M > 10~ 10 M Q ), a Pt < 10~ 20 ss^ 1 threshold 
of sensitivity should be reached. 



4. Solar-induced gravitational time delay 

The detection of dark planets close to Solar System is expected at low ecliptic latitudes. 
The SPSP deviation of the known planets and the Sun must be calculated and subtracted 
from the timing of the observed pulsars. The SPSP contribute from the Sun does not 
suffer the Re annual modulation because, with except for lower-order eccentricity effects, 
the observer-deflector distance Dd does not vary. Figure shows the Shapiro time delay 
A tocrai. induced by the Sun on the timing of a pulsar at an average distance of 4 Kpc 



(Taylor et al. 1993). The four curves correspond to four different values of the ecliptic 
latitude of the pulsar. The maximum Shapiro time delay is found when the position 
vector of the pulsar, D s , grazes the surface of the Sun and decreases at high ecliptic 
latitudes. Figures || and ^ show the corresponding residual time derivatives of the period, 
P. 

The formulation used in section is valid for impact parameters b = u m i n Re <C Dd 
and it is suitable to evaluate the maximum SPSP effect of the Sun. Taking characteristic 
values of Dd « 1 AU and (3± fts 10~ 4 we find: 




R B ~ V 2r s D d « 2.97- 10 9 J — — £- cm, (20) 



Re u„ 



9.91-10* = -£l \ Oh) 7^7 (21) 



Fo = w 9 - 94 ' 10 hod U^JUeJ h^J ss ss ' (22) 

The SPSP effect is huge (\P\ ~ 4.2-10~ 10 ss" 1 ) when 6 = u min R E = Rq {u min = 23.41). 
The characteristic time of the maximum SPSP effect is t c ~ 0.27 day. 
In order to avoid any additional and confusing refractive index delay (due to solar plasma 
even at large impact parameters) we must compare the two effects by comparing the per- 
turbations, Arig rav and An p , induced on the refraction index. The gravitational refraction 
index is: 

n grav ~l-^U=l + 4.22-lO- 6 (^ r )(^-) (23) 
while the plasma refraction index: 



M Q J \ R 
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where v v = J (e 2 n) / (Tim) = 8.98 • lO 3 ^/ n e Hz. For electrons, we have 



Aflg raV ^ f V \ ^ / Tip \ ^ 



r 



-l 



An p - im \lGHz) VlOScm-aJ \R Q ) ' (25) 
The density of electrons decreases with the distance from the solar surface (see figure 
|l0| ). At a distance of one solar radius from the surface of the Sun (r = 2R@, u m i n ~ 
46.8) n e ~ 10 5 cm,- 3 and therefore An grav /An p ~ 0.53 [v/ (1 GHz)] 2 . Close to Earth 
n e ~ 2 • 10 3 cm -3 and r = 1 and, again, An grav / An p ~ 0.49 [i// (1 GTJz)] 2 . From 
these values we deduce that, near Sun, pulsars might be better observed at high radio 
frequencies (y > lAGHz), which include the few known optical (X, 7) pulsars. In this 
high frequency range the SPSP effect of the Sun is dominant. 

5. Jupiter-induced gravitational time delay 

Jupiter is the most massive planet in our Solar System that can produce SPSP. Its 
geocentric orbit is a representative example of the annual Re modulation. Figure ^ shows 
the geocentric orbit of Jupiter during a twelve-year revolution. Because the distance of 
Jupiter from the Sun is about 5.2 AU, the geocentric distance Dd varies periodically from 
4.2 AU to 6.2 AU (A/Dd — 0.2), nearly twelve times during a revolution of Jupiter. 
The evaluation of the maximum SPSP effect Jupiter can produce can be made following 
the same procedure used for the case of the Sun. We find: 




/,V- V 2r,/), .-.2.11.1..^ (t— rrH . .- ,„,. (261 

t c = ^= « 7.15 • 10 3 ^~ )(-£-] s , (27) 

The equatorial radius of Jupiter is Rj = 7.14 • 10 9 cm, its mass is Mj ~ 0.78 • 10~ 3 Af Q 
and the average speed of revolution (3± ~ 4.36 • 10 -5 . Therefore, the average Einstein 
radius of Jupiter is Re — 1-89 • 10 s cm (u m i n — 37.82) and the maximum SPSP reaches 
a level of \P\ ~ 4.47 • 10 -14 s s _1 for a characteristic time of t c ~ 173 s. 
We also performed a simulation of the realistic SPSP produced by Jupiter on a distant 
pulsar. The pulsar is located at an average distance of 4 Kpc in the galactic centre direc- 
tion. Figure [ll] shows the gravitational time delay during a twelve- year cycle of revolution 
of Jupiter around the Sun, as seen from a geocentric observer. We can distinguish twelve 
peaks of gravitational time delay the magnitude of which increases as the position of 
Jupiter tends to the galactic centre direction. We can recognize the classic bell-shaped 
curve of a SPSP with a one-year modulation superimposed. Figure |l2| shows the corre- 
sponding residual time derivative of the period. The effect is not the maximum possible 
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because the minimum impact parameter of Jupiter with respect to the source is large 
compared with the Einstein radius (a little less than 1 AU). 



6. Mars-induced gravitational time delay 

When we take planets with orbits closer to the Sun then that of Jupiter, the period 
of revolution of gets comparable to one year and the modulation gets more strong and 
complicated in its time evolution. The case of Mars is shown in figure [l3|. Again, the pulsar 
is located in the galactic centre direction. The Einstein radius of Mars is much smaller 
than that of Jupiter while the minimum impact parameter reached in the simulation is 
practically the same. This implies a much larger u m i n and, so, a weaker event. 
Again, the maximum effect can be deduced from the equations: 

A.^.^u^)-^) (_-_)' (31) 

The mass of Mars is M M ~ 2.63- 1CT 7 M Q , the radius R M = 3.39- 10 8 cm (u mm = 294.78) 
and the average speed of revolution /3j_ ~ 8.05 • 10~ 5 . The maximum SPSP reaches a 
level of \P\ ~ 1.14 • 10~ 14 s s~ x for a characteristic time of t c ~ 228 s. 



7. Effects due to Earth rotation and Moon revolution 

We now check what kind of gravitational disturb both the moon and the rotation of 
the Earth should produce on the period time derivative of the pulsar. In the case of 
the Moon we can easily evaluate the maximum effect of the Shapiro effect from our 
usual formulas, when the trajectory of the signal grazes the surface of the moon. Setting 
Dd = 3.84 • 10 10 cm as the mean Moon-Earth distance , R m0 on = 1-74 • 10 s cm as the 
radius of the Moon and rs moon = rg^/81.3 as the Schwarzschild radius of the Moon the 
corresponding Einstein radius is RE moon — 2.89 • 10 cm. Taking R m0 on as the impact 
parameter u rn i„ RE moon we find 

P 0moon = 2rS(B = 4.29- 10 - 16 ss- 1 , (32) 

Rmoon 

which is quite a small effect. 

In addition we should take into account the gravitational perturbation produced by the 
rotation of the Earth. Let us take a simple limiting case. The pulsar is visible by the 
observer during the whole rotation and, at some instant of time, the trajectory of the 
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signal is orthogonal to the Earth radius (see figure |l4|) . In this simple case the maximum 
time delay during a rotation of Earth is: 



At rot ~ — £■ In 
c 



sin(2 ijj) + \Jl + sin(2V) 2 



(33) 



where ^ is the latitude of the observer. Suppose it be 2 7r/3, then At rot ~ 0.78 • 
r S e /c = 2.31 • 10 ~ n s and the mean residual time period derivative would be 
(P) ~ At rot / (0.5 doy) = 5.35 -lO" 16 sr 1 . 

8. Conclusions 

The Shapiro Phase might be able to discover not only dark matter but, once cali- 
brated, it may even lead to new mini planets discovers on our ecliptic plane. Their 
presence might be of relevance by tidal disturbances during rare encounters on Earth 



past and future history and life evolution (Fargion & Dar 1998). Finally recent EGRET 



data ( pixon ct al. 1998 ) on diffused Gev galactic-Halo might find an answer, among the 



others ( Fargion et al. 1998 ), by (~ pc) hydrogen molecular clouds interacting by proton 



(Gev) cosmic rays ( Dc Paolis et al. 1995| ). Microlenses at those large radii, (by Kirchoff 



theorem) are inefficient. Therefore the Shapiro Phase Delay might be the unique probe 
to verify such an evanescent dark matter barionic candidature. 

When taking in account the existence of celestial objects trapped by our Solar System we 
can view a one-year modulation on the SPSP produced by them as an outstanding unique 
feature. Being produced by celestial objects laying only near the geocentric observer this 
modulation can be used as a discrimination method. Among the possible methods of en- 
hancing the identification of some low-mass celestial objects within our Solar System we 
considered the observation of more than one pulsar at a time within a subset of pulsars 
near the ecliptic. A random combined statistic on many pulsars would deteriorate the 
predicted one-year modulation, because different modulations on many pulsars would 
interfere destructively. On the contrary a phase matching of the modulation may greatly 
increase by coherent interference the sensitivity, leading to the discover of the vectorial 
motion of the "dark" object. In any case, as a first point of view on this matter, we can 
say that the combination of both the methods could be taken into act. A one-year mod- 
ulated warning system would be useful for low-mass objects which can affect the residual 
timing of a pulsar at intensities well within the present observational detectability, while 
a combined statistic on many pulsars would be useful for the discovery of very low-mass 
celestial objects which would produce an SPSP of very low intensity. In view of a possible 
need of the discover of dangerous impact trajectories of asteroids, and in view of a need 
of monitoring incoming mini-planets, the cross-time correlation of the periods of a net of 
pulsars may offer a unique tool able to see in the dark. 
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Table 1. Pulsars with negative P. Comparison between the P in the two catalogues: 



(Taylor et al. 1993) fourth and fifth columns and (Taylor et al. 1995), sixth and seventh 
columns. 









P in the old catalogue 


P in the new 


catalogue 


pulsar 


Environment 


P (ms) 


P ( ss" 1 ) 


Epoch (MJD) 


P {ss- 1 ) 


Epoch (MJD) 


B0021-72C 


C 


5.8 


(-4 ±5) KT 17 


47858.521 


(-4.98 ±0.01) 10~ 20 


47858.5 


B0021-72D 


C 


5.4 




47858.521 


(-2.8 ±0.2) 10" 21 


48040.7 


B 1744- 24 A 


B , C 


11.6 


(-1.9 ±0.2) 10~ 20 


48270.0 


(-1.9 ±0.2) 10~ 20 


48270.0 


B1813-26 




592.9 


(-3 ±3) 10~ 16 


42004.1 


(±6.65 ±0.07) 10~ 17 


48382.0 


B2127+11D 


C 


4.8 


(-1.075 ±0.012) 10~ 17 


47632.52 


(-1.075 ± 0.012) 10" 17 47632.52 


B2127+11A 


C 


110.6 


(-2.107 ± .003) 10" 17 


47632.52 


(-2.107 ± .003) 10~ 17 


47632.52 



In the second right column the notations mean : C = globular cluster , B = binary pulsar 
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Table 2. Table of pulsars laying within an angular radius of <~ 6.38° close to the ecliptic. 



Pulsar B 


Pulsar J 


1 (deg) 


b (deg) 


P(s) 


P (ss- 1 ) 






0254-53 


0255-5304 


-90 


13 


-55 


31 





448 


(3.06 ± 0.03) 


10" 


17 


0447-12 


0450-1248 


-148 


92 


-32 


62 





438 


(1.033 + 0.008) 


10" 


16 


0450-18 


0452-1759 


-142 


92 


-34 


08 





549 


(5.7564 ± 0.0004) 


10" 


15 




0459-0210 


-158 


60 


-25 


75 


1 


133 


(0.00 + 0.00) 






0523+11 


0525+1115 


-167 


30 


-13 


24 





354 


(7.362 ± 0.008) 


10" 


17 


0525+21 


0528+2200 


-176 


14 


-6 


89 


3 


746 


(4.00321 ± 0.00006) 


10" 


14 




0533+04 


-159 


91 


-15 


32 





963 


(1.50 + 0.40) 


10" 


16 


0531+21 


0534+2200 


-175 


44 


-5 


78 





033 


(4.209599 ± 0.000002) 


10" 


13 




0538+2817 


179 


71 


-1 


68 





143 


(3.665 ± 0.005) 


10" 


15 


0540+23 


0543+2329 


-175 


63 


-3 


31 





246 


(1.542378 ± 0.000006) 


10" 


14 


0609+37 


0612+3721 


175 


45 


9 


09 





298 


(5.89 + 0.01) 


10" 


17 


0611+22 


0614+2229 


-171 


20 


2 


39 





335 


(5.9630 ± 0.0005) 


10" 


14 


0626+24 


0629+2415 


-171 


18 


6 


22 





477 


(1.99705 ± 0.00019) 


10" 


15 


0917+63 


0921+6254 


151 


43 


40 


72 


1 


568 


(3.6081 ± 0.0004) 


10" 


15 




1518+4904 


80 


80 


54 


28 





041 


(0.00 + 0.00) 








1627+1419 


30 


02 


38 


31 





491 


(3.93 + 0.03) 


10" 


Hi 


1633+24 


1635+2418 


42 


99 


39 


88 





491 


(1.194 + 0.002) 


10" 


Hi 




1640+2224 


41 


05 


38 


27 





003 


(2.9 + 0.2) 


10" 


21 




1645+1012 


27 


71 


32 


54 





411 


(8.13 + 0.02) 


10" 


17 


1709-15 


1711-1509 


7 


41 


14 


01 





869 


(1.104 ± 0.001) 


10" 


15 




1713+0747 


28 


75 


25 


22 





005 


(8.52 ± 0.02) 


10" 


21 


1718-02 


1720-0212 


20 


13 


18 


93 





478 


(8.7 + 0.2) 


10" 


17 


1717-16 


1720-1633 


7 


37 


11 


53 


1 


566 


(5.82 ± 0.02) 


10" 


15 


1718-19 


1721-1936 


4 


86 


9 


74 


1 


004 


(1.59 + 0.02) 


10" 


15 


1726-00 


1728-0007 


23 


02 


18 


28 





386 


(1.124 + 0.004) 


10" 


15 




1730-2304 


3 


13 


6 


02 





008 


(1.9 + 1.9) 


10" 


20 


1730-22 


1733-2228 


4 


02 


5 


74 





872 


(4.21 ± 0.09) 


10" 


17 


1732-02 


1734-0212 


21 


90 


15 


92 





839 


(4.2 + 0.2) 


10" 


Hi 


1732-07 


1735-0724 


17 


27 


13 


28 





419 


(1.21496 ± 0.00009) 


10" 


15 


1738-08 


1741-0840 


16 


95 


11 


30 


2 


043 


(2.274 ± 0.004) 


10" 


15 


1740-03 


1743-0337 


21 


67 


13 


41 





445 


(3.17 + 0.14) 


10" 


15 


1740-13 


1743-1351 


12 


69 


8 


20 





405 


(4.81 ± 0.06) 


10" 


Hi 




1744-2334 


4 


47 


2 


97 


1 


683 


(0.00 ± 0.00) 






1742-30 


1745-3040 


-1 


44 


-0 


96 





367 


(1.06649 ± 0.00001) 


10" 


14 


1745-12 


1748-1300 


14 


01 


7 


65 





394 


(1.212 + 0.002) 


10" 


15 


1745-20 


1748-2021 


7 


72 


3 


80 





289 


(4.00 ± 0.01) 


10" 


Hi 


1744-24A 


1748-2446A 


3 


83 


1 


69 





012 


(-1.9 + 0.2) 


10" 


20 


1744-24B 


1748-2446B 


3 


84 


1 


69 





443 


(0.00 ± 0.00) 
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Table 2. Continued 



1746-30 


1749-3002 


-0 


54 


-1 


24 





610 


(7.871 ± 0.002) 


10" 


15 


1747-31 


1750-3157 


-2 


01 


-2 


51 





910 


(1.97 ± 0.02) 


10" 


16 


1749-28 


1752-2806 


1 


53 


-0 


96 





563 


(8.1394 ± 0.0008) 


10" 


15 


1750-24 


1753-2502 


4 


25 





50 





528 


(1.413 ± 0.002) 


10" 


14 


1753-24 


1756-2435 


5 


02 


o 


04 


o 


671 


(2.8500 ± 0.0001) 


10" 


14 


1754-24 


1757-2421 


5 


30 





01 





234 


(1.30 ± 0.03) 


10" 


14 


1756-22 


1759-2205 


7 


46 





80 





461 


(1.08 ± 0.14) 


10" 


14 




1759-2922 


1 


19 


-2 


87 





574 


(0.00 ± 0.00) 






1757-23 


1800-2343 


6 


13 


-0 


12 


1 


031 


(0.00 ± 0.00) 






1758-23 


1801-2306 


6 


81 


-0 


07 





416 


(1.12983 ± 0.00003) 


10" 


13 


1757-24 


1801-2451 


5 


25 


-0 


88 





125 


(1.2790 ± 0.0001) 


10" 


13 


1758-29 


1801-2920 


1 


43 


-3 


24 


1 


082 


(3.292 ± 0.004) 


10" 


15 


1800-21 


1803-2137 


8 


39 


o 


14 





134 


(1.34327 ± 0.00001) 


10" 


13 


1800-27 


1803-2712 


3 


49 


-2 


53 





334 


(1.73 ± 0.03) 


10" 


17 




1804-2718 


3 


54 


-2 


82 





009 


(0.00 ± 0.00) 






1804-27 


1807-2715 


3 


84 


-3 


26 





828 


(1.225 ± 0.006) 


10" 


14 


1805-20 


1808-2057 


9 


45 


-0 


39 





918 


(1.7099 ± 0.0009) 


10" 


14 


1806-21 


1809-2109 


9 


41 


-0 


71 





702 


(3.822 ± 0.001) 


10" 


15 




1809-3547 


-3 


41 


-7 


84 





860 


(0.00 ± 0.00) 






1809-173 


1812-1718 


13 


10 





53 


1 


205 


(1.9119 ± 0.0004) 


10" 


14 


1809-176 


1812-1733 


12 


90 





38 





538 


(9.72 ± 0.08) 


10" 


16 


1813-17 


1816-1729 


13 


43 


-0 


42 





782 


(7.254 ± 0.006) 


10" 


15 


1813-26 


1816-2649 


5 


21 


-4 


90 


o 


593 


(6.65 ± 0.07) 


10" 


17 


1814-23 


1817-2312 


8 


48 


-3 


28 





626 


(0.00 ± 0.00) 






1813-36 


1817-3618 


-3 


19 


-9 


37 


o 


387 


(2.0465 ± 0.0014) 


10" 


15 


1817-18 


1820-1818 


13 


20 


-1 


72 





310 


(9.41 ± 0.03) 


10" 


17 


1819-22 


1822-2256 


9 


34 


-4 


37 


1 


874 


(1.353 ± 0.003) 


10" 


15 




1822-4210 


-8 


12 


-12 


85 





457 


(0.00 ± 0.00) 






1820-30A 


1823-3021A 


2 


78 


-7 


91 





005 


(3.385 ± 0.001) 


10" 


18 


1820-30B 


1823-3021B 


2 


78 


-7 


91 





379 


(3.15 ±0.03) 


10" 


17 


1820-31 


1823-3106 


2 


12 


-8 


27 





284 


(2.92294 ± 0.00018) 


10" 


15 


1821-19 


1824-1945 


12 


27 


-3 


10 





189 


(5.2248 ± 0.0013) 


10" 


15 


1821-24 


1824-2452 


7 


79 


-5 


57 





003 


(1.61845 ± 0.00005) 


10" 


18 


2151-56 


2155-5641 


-22 


95 


-47 


05 


1 


374 


(4.23 ±0.10) 


10" 


15 


2321-61 


2324-6054 


-39 


57 


-53 


17 


2 


348 


(2.60 ±0.11) 


10" 


15 
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x , image plane 




MACHO 



D ds = (1 - x d ) D s D d = x d D, 



Fig. 1. Relative position of source, MACHO , observer and primary image. As u — > oo, 
x + — > and the image corresponds to the source. v±_ accounts for the relative pulsar- 
MACHO-observer orthogonal motion velocity. 



1 




Umin * 10 


0.5 


^ ^^^^^^ 











0.5 






-1 







-20 -10 10 20 

ttc 



Fig. 2. F(t) as a function of t / t c . 
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PSR 




^ = ^ e 0- e = o 

D d D d 



Fig. 4. Relative position of the pulsar, Sun and Earth, ip is the ecliptic longitude and 
ip = 0, while 7r/2 — the ecliptic latitude. 




Fig. 5. Geocentric orbit of the Sun, Jupiter and the projection of the galactic plane. 
The line labeled as 7 indicates the vernal point direction, while the labels GC and GAC 
indicate the directions toward the galactic centre and anti-centre respectively. The orbital 
trajectories are calculated for a period of twelve years (UT) starting from 1998 January 
1. 
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Fig. 6. Mass of the deflector required to produce a P^-level deviation on the P of the 
pulsar. The angular impact parameter is the angular width in latitude of the subsample 
of pulsars shown in figure || and table H. The plot is a function of the distance Dd of the 
hypothetical planet from the center of the Solar System. 
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0.32 




1 2 3 4 5 6 

tp , rad 

Fig. 7. Gravitational time delay produced by the Sun on the timing of a pulsar laying 
at a mean distance ds = iKpc, plotted for four values of the orientation of the ecliptic. 
When = Tj- — •TCg"" the position vector of the pulsar, D 3 , grazes the surface of the Sun. 
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io- 12 
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Fig. 8. Residual period time derivative on the timing of a pulsar as a result of the 
gravitational field of the Sun for the same ecliptic latitudes and D s chosen for the plot 
shown in figure The complete plot for the maximum P deviation is shown in figure 0. 
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3.075 3.1 3.125 3.15 3.175 3.2 

(p , rad 



Fig. 9. Total plot of the maximum residual period time derivative as a result of the 
gravitational field of the Sun for the case $ = ? — R ^ un . 
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10 7 10 s 10 9 10 10 10 11 10 12 10 13 
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.0001 .001 .01 .1 1 10 100 



SOLAR RADII 



Fig. 10. Electron density, n e , from the solar surface. 
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0.00051 
0.0005 




0.00046 



2000 
d , day 



Fig. 11. Time delay as a result of the gravitational effect of Jupiter on the signal of an 
hypothetical pulsar laying at a distance of 4 Kpc in the galactic centre direction. 
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Fig. 12. Residual period time derivative for the plot in figure 
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200 400 600 800 1000 




200 400 600 800 1000 

d , day 



Fig. 13. Residual time derivative for a pulsar laying in the galactic centre direction at a 
distance of 4 Kpc, affected by the gravitational field of Mars, as observed by a geocentric 
observer (upper figure). The lower figure shows the same event with the modulation 
removed. The simulation covers the same twelve-year time interval taken in the case of 
Jupiter. 
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PSR 




Fig. 14. Path variation as a result of the rotation of Earth. 



